Bound of dissipation on a plane Couette dynamo 
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Variational turbulence is among the few approaches providing rigorous results in turbulence. In 
addition, it addresses a question of direct practical interest, namely the rate of energy dissipation. 
Unfortunately, only an upper bound is obtained as a larger functional space than the space of 
solutions to the Navier-Stokes equations is searched. Yet, in some cases, this upper bound is in 
good agreement with experimental results in terms of order of magnitude and power law of the 
imposed Reynolds number. In this paper, the variational approach to turbulence is extended to 
the case of dynamo action and an upper bound is obtained for the global dissipation rate (viscous 
and Ohmic). A simple plane Couette flow is investigated. For low magnetic Prandtl number Pm 
fluids, the upper bound of energy dissipation is that of classical turbulence {i.e. proportional to the 
cubic power of the shear velocity) for magnetic Reynolds numbers below and follows a steeper 
evolution for magnetic Reynolds numbers above P^^ {i.e. proportional to the shear velocity to the 
power four) in the case of electrically insulating walls. However, the effect of wall conductance is 
crucial: for a given value of wall conductance, there is a value for the magnetic Reynolds number 
above which energy dissipation cannot be bounded. This limiting magnetic Reynolds number is 
inversely proportional to the square root of the conductance of the wall. Implications in terms of 
energy dissipation in experimental and natural dynamos are discussed. 

PACS numbers: 47.27.N-, 52.30.Cv, 91.25.Cw 

Keywords: dynamo; energy dissipation; variational turbulence 
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I. INTRODUCTION 

Natural dynamos exist whenever the conditions of their 
existence is possible, i.e. when a sufficiently large mag- 
netic Reynolds is reached. One can imagine that this is 
an additional route for mechanical energy dissipation and 
that it is more likely for natural systems to take it than 
not to take it. It has been argued sometimes that a state 
of maximal dissipation rate should be reached and this 
idea has been used as a closure assumption for turbulence 
{e.g. Malkus Q). It is not necessary to make such an 
assumption, as explained clearly by Howard 0, and yet 
one can draw useful information from the determination 
of rigorous lower and upper bounds on energy dissipa- 
tion in turbulent flows. These bounds are obtained in a 
larger functional space than the solutions to the Navier- 
Stokes equation, hence they are not necessarily attained. 
However, in a number of cases, turbulent flows lead to 
a dissipation rate of the same magnitude as this upper 
bound m. 

The variational approach to turbulence was introduced 
by Malkus Q, Busse and Howard @. More recently, 
this approach was reformulated by Doering and Con- 
stantin [7| and expressed in a simpler way, using the con- 
cept of a background function (not necessarily the mean 
flow of turbulence) following Hopf Q. The objective is 
to bound energy dissipation under the constraint of hor- 
izontally averaged energy balance (for statistically plane 
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invariant configurations). In a series of papers, the ap- 
proach was improved by optimizing the spectral Lagrange 
parameter Q and the background function [H . The final 
bound is better than the nearly rigorous bound of Busse. 
Our objective here is not to exhaust these possibilities of 
optimization but rather to apply the general principle to 
a new configuration, i.e. a dynamo problem. 

The idea of applying a variational approach to a mag- 
netohydrodynamic flow has been applied already to a 
Couette and Poiscuille flow subjected to an applied trans- 
verse magnetic field [loj . This paper has been a source of 
motivation for the present work with two significant vari- 
ations. First, there is now no imposed magnetic field and 
secondly, magnetic boundary conditions are different. In 
the paper by Alexakis et al, magnetic disturbances are 
constrained to vanish at the boundaries, which does not 
correspond to a physically plausible situation. It is as- 
sumed here that the fluid domain is bounded by an infi- 
nite domain of electrically insulating medium (or with a 
conducting solid wall of finite thickness in between). 

A plane Couette flow configuration is considered and 
our objective is to find an upper bound to the total en- 
ergy dissipation when a prescribed velocity is applied. 
There is no applied magnetic field of external origin 
and the problem may look like a purely hydrodynamical 
one. However, the flow may support dynamo action (see 
[O, [13) the amount of dissipated energy must then 
take into account Joule dissipation. Upper bounds will 
be obtained as a function of two dimensionless param- 
eters, the Reynolds number and the magnetic Reynolds 
number. In addition, the effect of electrically conducting 
walls of finite thickness will be investieated. 
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Section |TT] provides details on the flow configuration, 
notation, dimensionless variables and equations. In sec- 
tion mil the principle of decomposition with background 
flow is presented. The horizontally averaged energy bal- 
ance is obtained in section IIVI and the expression for the 
total energy dissipation (Joule plus viscous) is given in 
section|Vl Energy dissipation bounds are obtained in sec- 
tion IVII and IVIII respectively, when velocity fluctuations 
and magnetic fluctuations are considered respectively in 
addition to the background flow. An improved bound 
is determined numerically in section IVIIII Section HXl is 
devoted to a discussion of the bounds obtained and their 
relevance to experimental and geophysical conflgurations 
and section |X] to the directions in which the variational 
approach could be extended to deal with more relevant 
models of dynamo. 

It may be useful to provide some guidance on how to 
read this paper. The Hopf-Doering-Constantin method is 
explicitly introduced and subsequent calculations of up- 
per bounds are also detailed explicitly. A reader with no 
prior knowledge of the method can check all results with 
pen and paper until the end of section IVIII Section IVIIII 
does not provide any fundamentally new result and can 
be ignored in a flrst reading. This section has required 
some standard numerical calculations of eigenvalues re- 
lated to the magnetic spectral constraint, and the method 
is only sketched. This section serves two purposes: first, 
it provides a better upper bound as it corresponds to a 
background function for which the spectral condition is 
only just satisfied (it is zero for a particular disturbance) 
and secondly, it provides a confirmation that the upper 
bound derived analytically in section IVIII is relevant as 
it obeys the same scaling law as the numerical bound at 
large magnetic Reynolds numbers. 
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FIG. 1: Plane Couette flow configuration 



electric current density j is simply the curl of the mag- 
netic field B. 

As shown on Fig. [1] the thickness of the fluid layer 
is 2H — 2E where E is the thickness of the electrically 
conducting layer of the wall. That thickness may be zero 
E — 0, in which case the wall is electrically insulating. 
In all cases, the overall thickness of the electrically con- 
ducting domain is 2H. For simplicity, the electrical con- 
ductivity of the solid conducting layer is identical to that 
of the fluid. Hence, at z = H ~ E and z = ~H + E, 
there is continuity of all three components of the mag- 
netic field. Moreover, the induction equation is identical 
in the fluid and in the conducting layer. Therefore, there 
is no boundary condition to consider at z = H — E or 
z = —H + E. The magnetic boundary condition, at 
z — ±H, is that the magnetic field in the electrically 
conducting domain is matched continuously to a curl-free 
and divergence-free magnetic field outside (see appendix 

S). 



II. PLANE COUETTE FLOW 
CONFIGURATION 

The dimensionless Navicr-Stokcs and induction equa- 
tions can be written 

+u-Vu = -Vp+2 xB + Rer^V^u, (1) 

ot 

— + u • VB = B • Vu -f i?„^V2B. (2) 

In the equations above, the dimensional length and veloc- 
ity scales are chosen to be half the distance, H, and dif- 
ferential velocity, U, between the plates (see Fig. [1]) while 
the quantities H/U, pU^, Uy/pll, y/pJ^iU/H are taken 
as dimensional scales for time, pressure, magnetic field 
and electric current density respectively. The dimension- 
less parameters are the Reynolds number Re = UH/v 
and the magnetic Reynolds number Rm = puU H. The 
symbols p, fj., v and a denote density, vacuum magnetic 
permeability, kinematic viscosity and electrical conduc- 
tivity respectively. In this dimensionless formulation, the 



III. BACKGROUND VELOCITY 
DECOMPOSITION 

A frame of reference {x,y, z) is chosen with x aligned 
with the direction of the imposed shear velocity and z 
in the direction perpendicular to the plates. Its origin is 
half-way between the plates. 

The velocity field u is written as the sum of a steady 
parallel profile ^{z)ex satisfying the boundary conditions 
$(±(1 — e)) = ±1 and of a field v with homogeneous 
boundary conditions v(a;, y, ±(1 — e), t) =0, i.e. 

u{x,y,z,t) = <P{z)e^+v{x,y,z,t) (3) 

It must be stressed here that $ - the so-called back- 
ground velocity profile - needs not be the average value 
of the full velocity field u, neither in the temporal nor in 
the horizontal space average sense. 

Substituting this decomposition into the Navier-Stokes 
and induction equations leads to evolution equations for 
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V and B for any particular choice of a background profile: 

— + v-Vv + $^ + Vz<^'ex = -Vp + j X B 

+i?e-i [$"e,. + V^v] , (4) 

dB dB 

— + V • VB + = B • Vv + S^^'e^r 
ot ox 

+R7n'y'B, (5) 

where $' and denote the first and second derivative 
of the background velocity profile <f> respectively. 



IV. ENERGY BALANCE OF FLUCTUATIONS 

The dot product of equations ([4]) and ([5]) with v and 
B respectively are integrated over V a large rectangular 
volume —L < x < L, —L < y < L and — 1 < z < 
1. They are then integrated over a long period of time 
T. The following space and time averaged equations are 
obtained: 

{<^'v,v,) = ((j X B) • v) - Re-' 

-Re-' (Vv : Vv) , (6) 
= ($'i3,i?,) + (j • (v X B)) - i?7„i (f) , (7) 

where the Poynting flux going out of the control volume 
has been assumed to be zero (there is no external source 
or sink, see appendix]^ and where the following global 
space-time average is introduced for any quantity /: 



< / >= lim lim 



1 



T ri rL rL 



T^oo L->oa L'^T 



J-Uo Jo 



f dxdydzdt. (8) 



The work of Lorentz forces ((j x B) • v) and electromo- 
tive work (j • (v X B)) arc equal and opposite, so that the 
sum of ([6]) and ^ leads to the following energy balance: 



{<S>'v,v,) + {^'B,.B,)-Re-'(<i> 



-Re-' (Vv : Vv) - R;,' (j^) = 0. (9) 



V. TOTAL ENERGY DISSIPATION 

Let us take pvU^ /H^ as a dimensional scale for energy 
dissipation per cubic metre. The average dimensionless 
energy dissipation due to viscous effects is (Vu : Vu) and 
the average dimensionless Joule dissipation is -Pm^(j^). 
When u is expressed using the decomposition ([3]), the 
total dissipation V = (Vu : Vu) -I- P^' (j^) becomes 

P = (Vv : Vv) + ($'2) + 2 /$'^\ + P-^' {f) . (10) 



Combining the energy constraint , so as to remove the 
linear term {^'dv^/dz), gives the following expression for 
the dissipation: 

V = - [(Vv : Vv) + 2Re {<i>'vxV,)] 

-P,-' - 2i?™ i'P'BxB,)] . (11) 

Dissipation is bounded by the background dissipation 
when both conditions (Vv : Vv) -1- 2Re {^'v^Vz) > 
and (j^) - 2R,n {^^'B^B^) > are satisfied for aU ad- 
missible vector fields v and B, i.e. divergence-free vector 
fields satisfying the appropriate boundary conditions. 
Those two conditions arc called spectral conditions in 
the framework of the background method because they 
can be treated as an eigenvalue problem and this is 
preceisely the method that will be followed in section 
IVIIII However, it is also possible to ensure that those 
conditions are satisfied using other methods and they 
will be dealt with using functional inequalities in the 
following sections IVII and IVIIl 



VI. VELOCITY SPECTRAL CONDITION 

Equation pT|) shows that the spectral conditions can 
be treated independently for velocity and magnetic dis- 
turbances. The velocity spectral condition consists in 
ensuring that (Vv : Vv) -I- 2Re {^'vxVz) remains positive 
for all velocity fields v. 

In order to satisfy the spectral constraint at large 
Reynolds numbers, it is convenient to choose a back- 
ground velocity profile which is uniform in most of the 
fluid, with a linear profile on each side so as to recover 
correct boundary conditions (see fig. [n|. The thickness 
S of these linear parts (virtual boundary layers) is a free 
parameter. The 'production term' 2Re {^'vxVz) is con- 
fined to these regions where it can hopefully be balanced 
by viscous 'dissipation' (Vv : Vv) for a sufficiently small 
thickness S. 

This problem is solved as follows (Docring and Gibbons 
[l^l, Alexakis et al. [l3l)- Using the condition v = at 
the lower wall, velocity at another position z in the fluid 
can be bounded as follows using Schwartz relationship: 
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(12) 



Combining with the corresponding equation for Vz and 
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using Young relationship leads to: 
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l-e 



Vv : Vvdz. 



(13) 



-l+e 



Integrating u^^w^^' over the lower boundary layer yields: 

-l+e+S ^ pl-e 

VxVy,<^'dz < - Vv : Vvdz (14) 



l+e 



l+e 



Averaging over the x and y direction and taking into 
account the upper boundary layer leads to: 



(15) 



The velocity spectral condition is then satisfied as soon 
as: 



S < Re-\ 



(16) 



In terms of energy dissipation ((TT|) . assuming magnetic 
effects to be absent, this background profile with suf- 
ficiently small boundary layers (|16p provides an upper 
bound: 



V < 2Re 



(17) 



VII. MAGNETIC SPECTRAL CONDITION 

In the previous section, it has been possible to find 
a condition on S fulfilling the spectral requirement for 
velocity perturbations. The corresponding task for mag- 
netic perturbations cannot be exactly similar as the prop- 
erty V = on the walls has to be changed into physically 
sound magnetic boundary conditions. We are consider- 
ing that the magnetic field inside the fluid and boundary 
domains should match a potential field outside decaying 
to zero at infinity (or possibly to a non-zero constant if a 
non-zero mean electric current density is allowed in the 
electrically conducting domain). 

We cannot easily find a simple expression for the 
boundary condition for or Bz in the physical space, 
however we obtain below a boundary condition for the 
{x,y) average product B^B^. This will then been used 
to bound {(^'BxBz) in terms of (j^). 

Let us denote the horizontal average over the (x, y) 
plane with an overbar: 



,f[z) = lim ^ 



"'0 



f{x,y,z)dxdy. (18) 



Outside the domain, the curl of the magnetic field is zero. 
So is the product of its y component with Bz averaged 
over X and y: 



Using integration by part, divergence-free condition for 
B, another integration by parts and the z component of 
V X B = leads to: 



Bz{dzB, - d^Bz) = dzB^Bz- B^dzBz, 
~ dzBxBz + B^dyBy, 



dzBxBz — BydyBx, 
dzBxBz + BydxBy, 
dzABl - 0, 



(20) 



As the product B^Bz is vanishing at infinite z, we con- 
clude that BxBz is zero at the boundary of the conduct- 
ing domain, i.e. when z = ±1. Strictly speaking, this is 
enough to see that it will be possible to bound (B^Bz) 
in terms of the mean square of the gradient of B using 
Poincare's theorem, and then in terms of the mean joule 
dissipation. A convenient way of obtaining such a bound 
is to repeat the calculations leading to (|^D|) while retain- 
ing a non-zero electric current density for —l<z<l. 
This leads to: 



dzBxBz = jyBz - jzBy = (j X B)„ (21) 

which is the expression of the x component of Lorentz 
forces in terms of the magnetic stress tensor, averaged 
on constant z planes. The product B^Bz can be evalu- 
ated inside the electrically conducting domain using the 
boundary condition B^Bz = at z = ±1. 

B^z{z) = J ^dzB^zdz', 

= (j X B)„ (22) 

which can be bounded as follows: 



BxBziz) < 



j2B2dz'. 




(23) 



As shown in appendix |Bl the x and y averaged square 
magnetic field in the electrically conducting domain is 
bounded as follows: 



B2(z)<2(f). 
Hence, equation (|23p can be written: 

\b;b:{z)\ < V2^iTi{f) 



(24) 



(25) 



BJdzB^-dM=0. 



(19) 



The magnetic boundary condition is rather loose com- 
pared to the velocity boundary condition which leads 
here to a bound proportional to the square-root of z -I- 1 
for the product B^Bz rather than a linear dependence 
for the velocity product VxVz. 

Using equation (l25ll it is now possible to bound 
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(B^B..^')\ < 



1 



B,B,{z)\dz 

l-e 



< 



< 



4^/2 



e+S 

y/udu. 



(26) 



The magnetic spectral condition (j^^ — 2R„^ {^'B,j.Bz) > 
is satisfied as soon as: 



8\/2i?„ 
36~ 



{e + 6) 



3/2 _ 3/2 



< 1. 



(27) 



This relationship is then used to express S is terms of vis- 
cous dissipation for the background flow, hence to deter- 
mine an upper bound for global dissipation. This bound 
is plotted on Fig. [2] for various values of the wall thick- 
ness e. The limiting case of vanishing wall thickness (or 
rather of an electrically insulating wall) can be easily de- 
rived from (j27p and corresponds to: 



8V2 Rm S^/^ 



< 1, 



(28) 



and since the dissipation of the background profile is 2/(5, 
we have: 



V < 



256 Rl 



(29) 
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FIG. 2: Bound on dissipation obtained from equations (|17|l for 
the velocity spectral condition ('plus' signs), and (|28|l for the 
magnetic spectral condition, for different wall conductances 
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One must take the higher upper bound derived from 
the velocity and magnetic spectral conditions in Fig. [2] as 
both conditions must be satisfied. It is straightforward 
from equations (|17p and (|27p that the magnetic spec- 
tral condition will take over the velocity spectral condi- 
tion when the magnetic Reynolds number is exceeding 
9/(128P„)~0.07P-i. 

It can also be derived from ((27|) that the critical mag- 
netic Reynolds number at which our upper bound esti- 
mate diverges towards infinity is Rm = l/(4'\/2e). 



VIII. NUMERICAL DETERMINATION OF 
UPPER BOUNDS 

The velocity and magnetic spectral conditions appear- 
ing in pT|) are tackled numerically in this section. The 
reason for this is not to achieve a better (i.e. lower) up- 
per bound, although this is the case. The main reason 
is to check that our analysis in sections IVII and IVIII pro- 
vides results that are of the correct magnitude and not 
gross overestimates of dissipation upper bounds. For a 
given background profile, the numerical analysis provides 
Reynolds and magnetic Reynolds numbers such that the 
spectral conditions are satisfied for all admissible velocity 
and magnetic fluctuations, and this cannot be improved 
as one particular velocity disturbance and one particular 
magnetic disturbance make the spectral conditions just 
zero. 

The method is sketched here: the spectral condition 
can be treated as a p roblem of energy stability for the 
background flow |13| . This is changed into an eigen- 
value problem and the game consists in finding a Reynold 
or magnetic Reynolds number such that the maximum 
eigenvalue is zero. The velocity and magnetic fluctua- 
tions are Fourier transformed in the x and y directions, 
while they are expanded using Chebyshev collocation 
polynomials in the z direction. The optimal magnetic 
disturbance is always found in the limit of zero wave num- 
bers kx and ky, whereas the optimal velocity disturbance 
is found for k^ = and ky finite and increasing with the 
Reynolds number (see Fig. (g])). 

For this numerical analysis, the background function 
is no longer piecewise linear. In order to avoid numerical 
difficulties due to discontinuities in the velocity gradient, 
a hyperbolic sine function class is chosen $ = sinh(z(5). 
In the limit of large dissipation, the background velocity 
is zero everywhere except in thin boundary layers where 
the velocity profile is an exponential function. Similarly 
to what has been done in sections IVTl and IVlIi the back- 
ground functions contain a single free parameter, i.e, the 
typical thickness 5 of the boundary layers. Due to the 
new profile of the background function, the associated 
dissipation takes a different form that can be derived an- 
alytically: 



1 - e + sinh (2i^) 



(30) 
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For a given choice of (5, hence dissipation, a lower bound 
of Reynolds and magnetic Reynolds number is sought. 
This is repeated for various values of 6 and the result- 
ing curve can be read as an upper bound of dissipation 
for each value of Reynolds or magnetic Reynolds number 
(see Fig. |3]and Fig. [3] respectively). Both spectral condi- 
tions can be plotted on the same figure once a magnetic 
Prandtl number is specified. 

It is clear that this numerical procedure leads to im- 
proved upper bounds - by a factor 10 approximately - 
compared to those obtained analytically, while similar 
trends are obtained. 




FIG. 3: Numerical magnetic upper bound for energy dissipa- 
tion for different wall thicknesses 
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FIG. 4: Numerical upper bound for energy dissipation from 
the velocity spectral condition and associated spanwise wave- 
number ky of marginal perturbation 



IX. DISCUSSION 

For the simple Couette flow configuration under in- 
vestigation, and for electrically insulating boundaries, it 
has been found that a classical 'hydrodynamical' upper 
bound for dissipation holds when Rm is smaller than 
around P^^'. the dimensional upper bound of dissipation 
per unit mass is of order U^/H. For larger magnetic 
Reynolds numbers though, another scaling is followed 
and the bound on dissipation becomes proportional to 
fPa'^h'U'^. This new bound is higher than the 'hydrody- 
namical' bound and is independent of the half-distance 
H between the plates. 

When the boundary plates have the same electrical 
conductivity as the fluid and a finite thickness e, the re- 
sulting upper bound is changed dramatically. The upper 
bound corresponding to the magnetic part of the spectral 
condition increases suddenly and diverges to infinity for a 
value of the magnetic Reynold number depending on the 
thickness of the walls. This critical magnetic Reynolds 
number scales as e~^/^. This value can be reached much 
before the change of regime between 'hydrodynamical' 
and (insulating) 'magnetic' branches of upper bound es- 
timates {R„i ~ Pm^)- 

It should not be too surprising that dissipation can be 
unbounded at some finite magnetic Reynolds number. 
Such a situation can be easily simulated with a BuUard 
dynamo model. This is a solid-rotation dynamo with no 
possible feedback of the magnetic field on the structure of 
the flow. Hence magnetic energy (and dissipation) grows 
without limit above the threshold of linear instability. As 
soon as electrical conductivity is restricted to the fluid 
(e = 0), this behaviour is no longer possible. For any 
finite magnetic Reynolds number, energy dissipation is 
certainly finite. Yet the scaling law for this dissipation is 
distinctly above the 'hydrodynamical' a priori estimate. 

Whether one is inclined to believe that upper bounds 
of dissipation are generally relevant to turbulence or not, 
there is one firm and strong result from the present work 
which is the extended validity of the Kolmogorov scaling 
up to Rm ~ Pm^- It might have been a reasonable guess 
to predict a distinct increase of the dissipation as soon 
as Rm exceeds a critical threshold of order unity. This 
analysis shows that this is not possible with electrically 
insulating boundaries. This certainly has implications on 
the design of future Couette-type dynamo experimental 
setups, depending on whether a small or large dissipation 
is sought. Although there is no mathematical statement 
on how dissipation is divided between mechanical and 
electrical phenomena, it might be a better choice to have 
electrically conducting boundaries, should one want to 
increase the ratio of Joule to viscous dissipation. 

It is appropriate to discuss here how energy dissipa- 
tion was affected by dynamo action in those three exper- 
imental setups where a self-sustained liquid dynamo was 
observed. They are the Riga (Latvia) [l4l. Karlsruhe 
(Germany) [31 and Cadarache (France) [l6| dynamos. 
From these references, a table can be shown (table IH), 
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TABLE I: Extra dissipation due to dynamo action slightly 
above threshold in experimental dynamos 





Riga [14J 


Karlsruhe \15} 


Cadarache \16^ 


7o above 
critical Rm 


6-7 


6.7 


30 


% extra 
dissipation 


10 


42 


15 - 20 



where the estimated extra dissipation when dynamo ac- 
tion is present is given (in percentage) for a certain mag- 
netic Reynolds number (expressed in percentage above 
the threshold) . the Karlsruhe dynamo experiment can be 
clearly distinguished by the large amount of extra dissipa- 
tion compared to the other two. This experiment is also 
one where differential velocity is forced deep inside the 
setup. For the Riga and Cadarache dynamos, a turbine 
or propellers are set in rotation near the outer boundary 
of the setup. There can only be a gross comparison to the 
pure Couette flow configuration presented here, but the 
case of thick electrically conducting boundaries is also 
putting shear velocities well inside the electrically con- 
ducting domain and results in a huge increase of the dis- 
sipation bound. Conversely, it cannot be ruled out that 
Riga and Cadarache dynamos are still obeying the same 
classical 'hydrodynamical' law of dissipation, similarly to 
the Couette flow with electrically insulating walls. The 
authors of the dynamo experiments themselves specify 
that the extra dissipation is obtained by difference of the 
actual dissipation and the extrapolated curve of dissipa- 
tion below threshold. Even though it has been used to 
estimate Ohmic dissipation by some authors, it should 
be stressed here that there is no theoretical reason why 
this extra dissipation should be the Ohmic contribution 
to dissipation. For instance, in the present work, we use 
the fact that there is an exact cancellation between the 
work of Lorentz forces and the electromotive work, but 
we have no access to either of these terms. 

We have concentrated mainly on the low magnetic 
Prandtl number case, because our primary interest lies 
in the Earth's core dynamics and in the relevant liquid 
metal experiments. Nevertheless the results we have ob- 
tained are equally valid for large magnetic Prandtl num- 
bers. It is however anticipated that our dissipation bound 
will be a gross overestimate at low values of the hydro- 
dynamic Reynolds number. In that case, the flow is the 
simple uniform shear flow between the moving walls and 
such a flow is not going to sustain dynamo action eas- 
ily. This perspective should not be ruled out as strong 
Lorentz forces might change the flow in such a way that 
it could drive a dynamo, however this transition - if it 
exists - will be severely subcritical. 

Saturation can somehow be addressed within this 
work. The total energy dissipation is an upper bound 
for Joule dissipation, which can be used through equation 
(IBll) to obtain an upper bound of the magnetic energy in 



a saturated regime of dynamo action. This bound cannot 
give realistic estimates near dynamo threshold, as dissi- 
pation is still dominated by the viscous contribution. At 
large magnetic Reynolds numbers, this bound is much 
larger than equipartition between kinetic and magnetic 
energy. 

In the derivation of velocity and magnetic spectral con- 
ditions, there are strong similarities. The only difference 
is due to boundary conditions for velocity disturbances 
and magnetic field disturbances. This mere difference in 
the boundary conditions leads to fundamentally different 
upper bounds for dissipation. Magnetic boundary condi- 
tions are sometimes not well treated in general textbooks 
(see, for instance, the energy stability of hydromagnetic 
flows in the otherwise excellent book by Joseph , vol- 
ume II). It is confirmed again in the present work that 
one should pay great attention to the relevant physical 
boundary conditions. 



X. PERSPECTIVES 

It is certainly useful to derive such upper bounds for 
dissipation in the field of magnetofluid dynamics, and 
more particularly regarding dynamos. When complex 
situations are considered (magnetic field, global rota- 
tion, ...) there is no reliable heuristic approach that 
can provide good estimates for dissipation. Ordinary hy- 
drodynamic turbulence can be roughly tackled by Kol- 
mogorov's theory but it is not safe to extend it to other 
types of turbulence. When the 'Doering-Constantin- 
Hopf ' background profile method can provided an upper 
bound, this is a solid reliable result. The bound can be 
rather constraining for instance in the case of a Couette 
flow with electrically insulating boundaries, at low Pm- 
it is not an obvious result that dissipation must remain 
similar to hydrodynamic dissipation up until ~ Pm^- 
There are actually a lot of other configurations for which 
the method can be applied. 

Upper bounds for dissipation can be calculated for ex- 
perimental setups involving liquid metal flows. In Greno- 
ble, LGIT laboratory, we have the DTS setup involv- 
ing 40 liters of sodium in an imposed magnetic field 
fisl . [20| . Upper bounds for dissipation would provide 
some information about the importance of size scaling of 
the experimental setup. This is crucial when planning a 
larger setup that might sustain dynamo action. 

Thermal convection can also be taken into account to 
derive dissipation bounds. This is the driving force for 
the flow of liquid iron in the Earth's core. Upper bounds 
of dissipation may prove useful in terms of the thermal 
budget of Earth throughout its history. Although the 
configuration of the core of the Earth is not similar to 
the simple Couette flow considered here, it could be 
very relevant to examine the role of a thin electrically 
conducting layer at the bottom of the mantle, as it could 
make a large change in the upper bound for dissipation. 
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APPENDIX A: MAGNETIC BOUNDARY 
CONDITIONS AND ASSOCIATED POYNTING 
FLUX 

Outside our electrically conducting domain —H<z< 
H (i.e. the fluid layer and electrially conducting part of 
the plates), it is assumed that there exist no other elec- 
trically conducting domains, and no regions with mag- 
netic properties: let us just think of it as empty space. 
Coherent with the classical magnetostatic approximation 
used in the electrically conducting domain to derive the 
induction equation the displacement current is also 
neglected outside, so that the magnetic field satisfies the 
following Maxwell equations: 



V • B = 



V X B = 0, 



(Al) 



from which it follows that the magnetic field is harmonic. 
The magnetic field can be written as a Fourier integral: 



showing that it must also decay exponentially. Hence, the 
Poynting flux density E x B//^ decays exponentially and 
is thus as close to zero as one wishes some distance away 
from the fluid layer. From the assumption of statistical 
stationarity, there can be no accumulation of magnetic 
energy, hence this Poynting flux E x B//i must vanish 
when averaged along x and y, at any position z out- 
side the electrically conducting domain. This reasoning 
is valid for each wavenumber and the total Poynting flux 
is the sum of individual wavenumber components, hence 
this proves that the averaged Poynting flux must be zero 
at any position z and in particular at the boundary z = H 
and z = —H (or z = ±1 in dimensionless coordinates). 
Let us insist that this result is true when the Poynting 
flux is averaged on a plane of constant z (not pointwise) 
and when it is also averaged in time, assuming stationary 
turbulence. 

It should be noted that the above reasoning fails in 
two circumstances. One is when there are other electri- 
cally conducting domains. In that case, magnetic energy 
can be dissipated or generated in each domain and there 
can be a net exchange of energy between the fluid layer 
and the other domains. The other possibility is that en- 
ergy is radiated away from the fluid layer and goes away 
infinitely far. This is only possible when the displace- 
ment current is taken into account in Ampere's equation 
(|A1[) . Then a significant fraction of energy can be radi- 
ated away only when the timescale r and lengthscale I of 
the magnetic field are such that 1/t is comparable to the 
speed of light [2l| . This is a necessary condition for our 
system to work as an antena, and we shall not consider 
this limiting case. 



B 



B(fc^,fcy,z,t) e(''=-^+''='''')dfc^ dky 



(A2) 



and the harmonicity of B implies that each component 
(kx and ky fixed) must satisfy: 



-k^B, 



(A3) 



APPENDIX B: BOUND ON MAGNETIC 
ENERGY FOR A GIVEN AMOUNT OF JOULE 
DISSIPATION 

We show here that magnetic energy in the conduct- 
ing domain is bounded pointwise by the integral Joule 
dissipation. In our dimensionless terms, the following 
relationship hods for any value of z: 



where k'^ ^ k'^ + ky . As a consequence B must decay ex- 
ponentially away from the fiuid layer on a typical distance 
l/k. There is one exception when there is a significant 
magnetic contribution near fc = 0. However, we shall ex- 
clude that possibility on physical grounds: in that case, 
the exponential decay length increases without limit and 
there would be an infinite amount of magnetic energy 
stored oustside the fluid. That would take an infinitely 
long time to settle and is thus incoherent with our as- 
sumption of the existence of a stationary flow solution. 

Faraday's equation (another equation from the set of 
Maxwell's equations) is used to determine the electric 
field: 



V X E -■ 



dB 



(A4) 



B2(z)<2(f) 



(Bl) 



The magnetic field B is decomposed in poloi'dal-toroi'dal 
contributions: 

B = V(7xe,, + Vx (Vp x e^) + boxiz)ex + boy{z)ey (B2) 

The assumption of spatial statistical invariancc along 
X and y directions allows us to write the poloi'dal and 
toroidal contributions of B in terms of a Fourier integral 

[p, q]^ 11 [P, Q] {kx,ky,z,t) e^'''^-+'''yy'>dkxdky 



(B3) 

Each Fourier component can be considered individually, 
as magnetic energy and ohmic dissipation will be iust 
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the sum of energy and dissipation of each contribution. 
It is enough to prove (|Bip for aU contributions. The 
contribution boxiz)ex + boy{z)ey will be dealt with at the 
end of this appendix. Poloi'dal and toroidal contributions 
lead to 



B 



-dyQ - dxzP 



(B4) 



j = V X B = 



-dxV^p-d^.q 



(B5) 



Hence, x and y averaged magnetic energy and ohmic dis- 
sipation can be written 



B^z) = ^\Q\' 



y|5.P|Vy|P|' (B6) 



F(z) = ^ \-k'P + d',P\' + ^ + ^ IQI^ (B7) 

Integration by parts leads to the following expression 
for the global dissipation using the condition that the 
poloi'dal and toroidal components must decay far away 
from the fluid layer 



0^) = 



1 

-1 



pdz 



OO 



'—\dlP\'' + k*\d,P\'' + — \P\^dz 



/OO 7,2 1,4 
^ylS.gi' + ylQl'd^ (B8) 

The toroidal scalar function Q must vanish at z = ±1. 
Hence: 



Q{z) = J d.Qdz 
Hence by Cauchy-Schwartz 

IQl' < j dzj \d,Q\^dz 



(B9) 



< 2 / \d,Q\ dz 



(BIO) 



The treatment of the poloi'dal is slightly more involved 
as the boundary conditions available are less straightfor- 
ward: dzP ± P = at 2 = ±1. They can be used as 
follows 

2kP = {kP - dzP) + {kP + dzP) 

= I kd:,P - dl,Pdz' - I kd^P + dl^Pdz' CQll) 



from which the modulus of P can be bounded using the 
triangle and Cauchy-Schwartz inequalities 

2k\P\ < f k\d,P\ + \dl,P\dz' 



< 



12 J k'^\d^P\^dz' 



Hence 



k^ ipr < k' 



\d.P\^dz' 



12 I \dl^P\Uz' (B12) 



dl^P\ dz' (B13) 



A similar treatment is made on dzP 

2d,P = - {kP - d,P) + [kP + d,P) 
leading to a similar result 



(B14) 



\d,Pf<k''j \d,Pfdz' + J \dl,P\^dz' (B15) 



All three terms in (jB6[) can be bounded using (jB10[) . 
(jB13[) and (jBlSp and then compared to the contributions 
of the ohmic dissipation in (jBSp , leading to the expected 
result: 



B2(z)<2(f) 



(B16) 



Importantly, this result is independent of fc, so that it 
will equally apply to a sum of contributions of different 
wavenumbers kx and ky, hence proving the expected re- 
sult for the fields B and j. 

The horizontally independent contributions bQx{z)ex + 
boy{z)ey must now be considered. Their associated ohmic 
dissipation can be written 





dbox 


2 

+ 


dboy 




dz 




dz 



dz 



(B17) 



If the global electrical current flux in the x and y direc- 
tion is not zero, then this is an unphysical situation as 
magnetic energy in the space above and below the fluid 
would be infinite, and such a situation would take an in- 
finitely long time to be established. If it is zero, then box 
and boy vanish on z = ±1, and one can derive a bound 

for |B|^ using Cauchy-Schwartz inequality 



\box\'^ (z) 



dbpx 
dz 



dz 



< 2 



db 



Ox 



dz 



dz (B18) 



A similar inequality holds for &oyj so comparing (|B17[) 
and (IB18I) is also compatible with the inequality to prove. 
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